我正在尝试在R中运行一个间隔审查的COX模型。CENREG软件包具有一个函数IC_SP来执行此操作。但是,当我将其作为预测因子的因素变量时,我会遇到问题。data=data.frame(R=c(2,1,4,7,4,6,8,8,12,3),L=c(3,1,4,7,5,6,8,9,12,3),gender=c("m","f","f","m","m","f","f","m","f","m"))s=survival::Surv(data$R,data$L,type='interval2')icenReg::ic_sp(s~relevel(gender,"f"),model='ph',bs_sample